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Abstract 

We have established that the mean photon number obtained employing the conventional 
Hamiltonian for the twin light beams with the same frequency, produced by a subharmonic 
generator, is half that of the twin light beams with different frequencies. In view of this, we have 
found it to be quite appropriate to describe the process of subharmonic generation by the same 
Hamiltonian regardless of whether the twin light beams have the same or different frequencies. 
Applying this Hamiltonian, we have calculated the photon statistics and quadrature squeezing 
for the signal-signal and the signal-idler modes. Obviously, this time the mean photon number 
of the signal-signal modes is exactly the same as that of the signal-idler modes. In addition, we 
have found that the maximum global quadrature squeezing of the signal-signal or the signal-idler 
modes is 50% below the vacuum-state level and the maximum local quadrature squeezing of the 
signal-signal modes is 74.9% below the same level. 
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1 Introduction 

In subharmonic generation a pump photon of frequency u is down converted into a pair of highly 
correlated photons of frequencies wi and u} 2 , with wi and u }2 being the same or different. The 
statistical and squeezing properties of the twin light beams with the same or different frequencies 
(each light beam consisting of one photon from each pair) have been investigated by several authors 
[1-4]. It is found that the twin light beams are in a squeezed state, with the maximum quadrature 
squeezing being 50% below the vacuum-state level. We recall that the twin light beams with the 
same frequency are represented in the conventional Hamiltonian by and dF. For a given pump 
mode, we anticipate the mean photon number of the twin light beams with the same or different 
frequencies to be the same. However, the mean photon number of the twin light beams with the 
same frequency obtained, in the Appendix using the conventional Hamiltonian, is found to be half 
that of the twin light beams with different frequencies, as determined in section 4. We maintain the 
standpoint that the number operator d^d holds only for a light mode represented in the pertinent 
Hamiltonian by first-order annihilation and creation operators. The unexpected result we have 
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obtained for the mean photon number of the twin light beams with the same frequency must 
then be due to representation of these twin light beams by second-order annihilation and creation 
operators. 

We may envisage subharmonic generation as a process in which a three-level atom absorbs a pump 
photon and makes a transition from the bottom to the top level. Then it emits a photon and decays 
to the intermediate level. Finally, the atom emits another photon and decays to the bottom level. 
We certainly expect several three-level atoms to be involved in this process. The photons emitted 
from the top level constitute one light mode and those emitted from the intermediate level form 
another one. The two light modes may have the same or different frequencies. It therefore appears 
to be quite appropriate to describe the process of subharmonic generation leading to the creation 
of twin light modes with the same or different frequencies by the Hamiltonian [5] 

H = —b) + iX(b^aia2 — (1) 

where ai and 0,2 are the annihilation operators for the light modes emitted from the top and 
intermediate levels, b is the annihilation operator for the pump mode, A is the coupling constant, 
and ^ is proportional to the amplitude of the coherent light deriving the pump mode. We may 
refer to a Hamiltonian of the form described by ([T|) as first-order Hamiltonian. We assume that the 
operators oi and 02 commute and satisfy the commutation relation 

[fii, all = [02,4] = 1- (2) 

We next seek to calculate, applying the Hamiltonian given by Eq. m, the mean and variance of the 
photon number, the global quadrature squeezing, and the photon number distribution for the twin 
light modes with the same or different frequencies. We also wish to obtain the local quadrature 
squeezing for the twin light modes with the same frequency. 


2 Operator dynamics 



pump mode 






light mode oi 


light mode 02 


Fig. 6.1 A subharmonic generator. 


We consider here the case in which the twin light modes and the pump mode are in a cavity coupled 
to a vacuum reservoir via a single-port mirror. It is not hard to realize that 


d 


— (ai) = -iTr{p[di,H]) -b yrr((2ai/5a| - d\dip - pd\di)di), 


^y(a 2 ) = -iTr[p[d 2 ,H]) -b ^Tr{( 2 d 2 pd\ - d\d 2 P - ^ 0 ^ 02 ) 02 ), 
y(a|ai) = -iTr[p[d\di,H]) + ^Tr[{2dipd\ - d\dip - pajai)4ai), 


(3) 

(4) 


( 5 ) 








^{a\a2) = -iTr[p[a\a2, H]) + ^Tr[{2a2pa\ - a\a2P - pa\a2)a\a2), 


( 6 ) 


^{aia2) = -iTr[p[aia2.,H]) + '^Tr(^ 2 aipa\ - a\aip - pa\ai)aia‘^ 
+^rr ( [2a2pa\ - a\a2P - pa\a2)aia2 j , 


(7) 


in which ki and K2 are the cavity damping constants for light modes oi and 02, respectively. Now 
in view of Eq. © and the fact that 


Yrr((2aipa| - a\aip - / 5 a|ai)ai) = — 
^Tr{{2a2pa\ - a\a2P - pa\a2)a2) = -^( 02 ), 


-^Tr{^ 2 aipa\ — a\aip — / 5 d|ai)a|oi) = —Ati(ajai), 


H2 


^ Tr{{2a2pa\ - a\a2P - pa\a2)a\a2) = -^ 2 ( 0202 ), 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


we readily get 


^Tr\^ 2 aipa\ — a\aip — pa\ai)aia2 

+^Tr(^2a2pa\ — a\a2P — pa\a2)aia2 
= + «^2)(aia2), 


^(«i) = -^Ki(ai) - \{hal), 


^(62) = -^^2(02) - A(6a{), 


^(o|di) = —Ki(a|di) — A(Sldia2) — A(6a|a2), 
dt 


(0^02) = -K 2 {dld 2 ) - A(6laia2) - A(6a|a^), 


^(0102) = + K2)(dia2) - A(6ajai) - A(6a^a2) - X{b). 

On taking ki = k,2 = k, the steady-state solutions of the above equations are found to be 

(ai) = - — {bdl), 


(02) = - — {bd\), 

Kj 


(ajai) = --(6laia2) - -(6a|d^), 

Kj Ki 

(0^02) = --(6laia2) — —{bd\dl), 

K K 

(0102) = — — (hajai) — —(60202) — —(b). 

K K K 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 


Upon dropping the noise operator and in the absence of subharmonic generation (A = 0), one can 
write the quantum Langevin equation for the operator b as 


db 1 ^ 


(23) 

where k is the cavity damping constant. The steady-state solution of this equation is 


6 = ^. 


(24) 

Now on introducing (1241) into Eqs. (I18l)-(l22p. we arrive at 



(a,) = -^(4>, 


(25) 

(as) = 

Ki 


(26) 

(a|ai) = --(aios) - -{d\dl), 

K K 


(27) 

(a^as) = --(aios) - -(a|4)> 

K K 


(28) 

(aias) = --(ajai) - -(4«2) - 


(29) 

in which e is defined by 

2^X 

K 


(30) 

Applying Eqs. (12^ and (l26l). we easily find 



(oi) = ( 02 ) = 0. 


(31) 

Moreover, using Eqs. (ITTl). (f25|l. and (f2U|l. one gets 



(ai«i) k2_4^2’ 


(32) 

(4«2) = {aldi), 


(33) 

and 

It can also be readily verified that 


(34) 

(“ 1 ) = (« 2 ) = ( 0102 ) = 0. 


(35) 

We take 

0 = Q-l -|- 02 


(36) 

to be the annihilation operator for the superposition of light modes ai 
subharmonic generator. One can easily check that 

and 02 , produced by the 

[ 0 , a'*'] = 2. 


(37) 

We realize that the superposition of the two light modes, with the same or different frequencies, 
constitutes a two-mode light. We wish to call the superposed light modes with the same frequency 
the signal-signal modes and the superposed light modes with differen frequencies the signal-idler 
modes. It also proves to be convenient to refer to the subharmonic generator which produces the 
signal-signal modes as the degenerate subharmonic generator and the one which produces the signal- 


idler modes as the nondegenerate subharmonic generator. Finally, we would like to mention that 
the results described by Eqs. are valid for the signal-signal or the signal-idler modes. 





3 The Q function 


We next seek to obtain the Q function for the two-mode subharmonic light (the signal-signal or the 
signal-idler modes). This Q function is expressible as 

Q{ai,a2,t) = ^ J (fz(fr]4>a{z,ri,t)exp{z*ai-zal + rj*a2-rja2), (38) 

in which the antinormally-ordered characteristic function 4>a{z,ri,t) is defined in the Heisenberg 
picture by 

=rr(p( 0 )e-^*“iWe^“^We-’?*“ 2 We»? 4 (i)). ( 39 ) 

Employing the identity 

^ ^A+b+^[A,b] ^ 40 ^ 

and taking into account the fact that ai and 0,2 are Gaussian variables with zero mean, Ea. (l39p can 
be put in the form 


1 1 

4 >aiz,r],t) = exp[--{z*z + r]*r])exp[-{{za\ - z*ai + q&l - 7^*02)^)]. ( 41 ) 

It then follows that 

(j)a{z,r],t) = exp[-^{z*z + r]*r])]exp[{ — ^z*z{a\ai + aia|) -|- ^z*^al 
+ \z‘^a^^ - \p*p{a\a 2 + 020^) -|- \ni*‘^a 2 + \'rf‘a ^2 
+z*vi*aia2 + zria\a2 — zr]*a\a2 — 2 * 770102 )]. (42) 

Now on account of (l32]l . ([33|) . (p^ . and (f35]l . this equation goes over into 

0a(2,77) = exp [-0(2*2 -I- ri*r]) - b{zr] + 2*77*)] , ( 43 ) 


where 


and 


_ _ 2^2 
“ k 2 _ 4g2 


(44) 


b 


K£ 

K? — 4^2 


(45) 


Finally, upon introducing (1431) into (1381) and carrying out the integration, the Q function for the 
two-mode subharmonic light is found to be 


<5(0:1,02) 


\u^ — t;^] exp [—n(oi oi -|- 0902) 

vr^ L j 


77(0102 “k <^10^2)1 5 


(46) 


in which 


a 

= ^23^’ 

b 

a? — 52 


(47) 

(48) 


4 Photon statistics 


In this section we wish to calculate the mean photon number, the variance of the photon number, 
and the photon number distribution for the signal-signal modes as well as the signal-idler modes. 






4.1 The mean and variance of the photon number 


We define the mean photon number of the two-mode subharmonic light by n = (a^a). Then using 
Eq. ([36]) and taking into account Eq. ([351) . we easily find 

n = {a\ai) + (0^02), ( 49 ) 

so that in view of ()32|) and (j33|) . there follows 

This represents the mean photon number of the signal-signal or the signal-idler modes. We observe 
that the mean photon number obtained in the Appendix is half of the result given by Eq. (|5UI) . 

On account of Eq. (|37p . the photon-number variance of the two-mode subharmonic light, defined 
by 

(An)^ = ((a^a)^) — (51) 

can be put in the form 

(An)^ = -|- 2n — n^. (52) 

Now applying the fact that a is a Gaussian variable with zero mean, we get 

(An)^ = 2n-|--|-(a^^)(o^) (53) 

and on taking into account (|36l) along with (13511 . we arrive at 

(An)^ = 2n-|--|-4(a|a2)(aia2). (54) 

, the photon-number variance of the two-mode subharmonic 


Hence in view of Eqs. (|50p and 
light takes the form 

(An)^ = 






166^ 




4:K^£‘ 


2^2 


2 _ 4^2 _ 4£2^2 ^ 1^2 _ 4£2^2 


(55) 


In addition, we note that the equation of evolution of the mean photon number for the pump mode 
can be written as 


^{h^b) = — H]) -|- ^KTr[{2b0 — Pbp — 

OjL Zi 


Then using Eq. © and the fact that 


we readily get 


dt 


^KTr[{2bpP — b^bp — pPb)b%] = —K{b^b), 


(b^b) = —K{Pb) + pb + pb^ + X{b^did2) + A(6a|a2). 


The steady-state solution of this equation is 


(Pb) = —b + —6'*' -I- -(6'*'aia2) -|- -{bd\dl), 


(56) 

(57) 

(58) 

(59) 


so that in view of Eqs. 


and (|34l) . the mean photon number of the pump mode takes the form 


- 4e^ 


(60) 


The first term represents the mean photon number of the pump mode in the absence of the subhar¬ 
monic generation and the second one represents the mean photon number of light mode oi or light 
mode 02 . This is exactly what we would expect the mean photon number of the pump mode to be. 







4.2 The photon number distribution 


We finally proceed to calculate the photon number distribution for the signal-signal or the signal- 
idler modes. The probability for observing m + n signal photons or the probability for observing m 
signal photons and n idler photons is expressible in the form 


P{m, n) = 


d2r, 


a2r, 


■ [Q{ai , a2)e°‘*°'^+°‘2°^'^]\ai=a*=a2=a*=0- 


(61) 


Thus with the aid of Eqs. ([l6l) and (f6T]l . one can write 


P{m, n) - 


\y? — 


Q2r, 


m\n\ daf'da'^ da^da^ 

X exp [(1 - u){a\ai -h a 2 « 2 ) - v{aia 2 -h ata 2 )]«i=aj=a 2 =«^=o > 

so that on expanding the exponential functions in power series, we have 

\v? — v‘^^ (—1)^+^(1 — 


mini 


ijki 




X 




1 


a. 


Qin 


dafda\^ dd^daf 

- ai= 0 j= 02 =“ 2 “e 


(62) 


(63) 


Now upon performing the differentiation and applying the condition ai = a\ = a 2 = 012 = 0, one 
gets 


P{m, n) 


y - u"] 

mini 

{i + k)i {i + e)i {j + k)i {j + e)i 
{i + k — m)l {i + l — m)l {j + k — n)l {j + £ — n)l 


E 

ijki 


(_l)fc-|-^(l _ uy+iy^+^ 


x6, 


i+k,m^i+£,m^j+k,n^j-\-i,n 


(64) 


We note that k 
the form 


i = m — i = n — j. Therefore, for m = n the photon number distribution takes 


P{n,n) = [u^ — ^ 

j=0 


n!^(l — 

min-j)iy 


(65) 


We observe that Eq. (I65p represents the probability to observe 2n signal photons or n signal photons 
and n idler photons. 


5 Quadrature squeezing 

In this section we seek to obtain the global and local quadrature squeezing for the two-mode sub¬ 
harmonic light. 


5.1 Global quadrature squeezing 

Here we wish to determine the global quadrature squeezing (in the entire frequency interval) for 
the two-mode subharmonic light. We recall that the variance of the plus and minus quadrature 
operators for a two-mode light is given by 

{Aa±f = (a±,a±). 


( 66 ) 












where 


a+ = + a, (67) 

a_ = i(a^ — a), (68) 

and a is the annihilation operator for the two-mode light. Now on account of (|31l) along with (j36l] . 
(1671) . and (l68l) . we see that 

(a±) = 0. (69) 

Hence in view of this, Eq. (|66jl is expressible as 

(Aa±)2 = (a|). (70) 

We note that 

(Aa±)^ = (d^d) -|- (aa^) ± (d^ -|- d^^^) (71) 

and applying (|37l) the quadrature variance can be put in the form 

{Aa±f = 2 + 2{a^a) ±{ 0 ."^+ a^^). (72) 

With the aid of ((Ml) , one can rewrite Eq. (|72]l as 

(Au^)^ = 2 -)- 2(d|di d2d2) i (d^ -1- d| -1- d^ -b d^ -b 2did2 -b ^d^d^). (73) 

Therefore on account of Eqs. (f32]l . ([33|) . ([3^ . and (f35]l . the quadrature variance takes the form 

and 

(Aa_)2 = 2 + ^^. (75) 

Kj ZiS 

We immediately note that the two-mode cavity subharmonic light is in a squeezed state and the 
squeezing occurs in the plus quadrature. In addition, we note that for k = 2e the variance of the 
minus quadrature diverges. We then identify k = 2e as the threshold condition. 

Upon setting e = 0 in Eqs. (1711) and (f75]l . we find 

(Aa±)2 = 2. (76) 


We thus see that for e = 0 the cavity light is in a two-mode vacuum state in which the uncertainties in 
the two quadratures are equal and satisfy the minimum uncertainty relation. We wish to calculate 
the quadrature squeezing of the two-mode cavity subharmonic light relative to the quadrature 
variance of the two-mode cavity vacuum state. We therefore define the quadrature squeezing of the 
two-mode cavity subharmonic light by 


o 2-(Aa+)^ 
2 

Now combination of ((Till and (|77)) yields 


S 


2e 

K, “b 2e 


(77) 


(78) 


This represents the global quadrature squeezing of the cavity signal-signal or signal-idler modes. 
We note that at steady state and at threshold there is a 50% quadrature squeezing below the 
vacuum-state level. 






We next seek to obtain the quadrature squeezing of the two-mode output subharmonic light. To 
this end, we demand that the quadrature variance of the two-mode output subharmonic light must 
be the sum of the quadrature variance of the transmitted two-mode cavity subharmonic light and 
the quadrature variance of the two-mode reflected input modes. One can then write 

= K{Aa±f -h (1 - K)(Aain±)^, (79) 

so that on account of (1711) and (f75]l and the fact that 

{Aain±)^ = 2, (80) 

we easily find 

(Aa°“*)2 = 2- (81) 

and 

A CT 

(Aa°“*)2 = 2 +-. (82) 

K -|- 2e 

Moreover, the quadrature variance of the two-mode output vacuum state can be written as 

(Aa±“*)^ = K{Aa±)l + {1 - K){Aain±f. (83) 

Hence in view of (I76p and (I80p . there follows 

(Aa^“*)2 = 2. (84) 


Now we define the quadrature squeezing of the two-mode output 


i^OUt _ 


2 - {Aa°^*f 


(85) 


Then with the aid of (j8ip . we obtain 


_ 


2Ke 


( 86 ) 

K + 2e ^ ^ 

Finally, at threshold and for k = 0.8, we see that the qnadrature squeezing of the two-mode output 
subharmonic light is 40% below the vacuum-state level. We thus note that the quadrature squeezing 
of the two-mode ontput subharmonic light is less than that of the two-mode cavity subharmonic 
light. This must be due to the vacuum reservoir quadrature fluctuations. 


5.2 Local quadrature squeezing 


We finally seek to calculate the local quadrature squeezing (in a given frequency interval) for the 
signal-signal modes. To this end, we define the spectrum of quadrature fluctuations for the cavity 
signal-signal modes with central frequency wq by 

S±{uj) = — {a±{t),a±{t + T))sse'^‘^~‘^°^'^dT. (87) 

J — OO 

Upon integrating both sides of (IHTP over u, we easily get 


S±{uj)du} = (Aa±)^, 


( 88 ) 


with (Aa±)^ being the steady-state global quadrature variance. On the basis of this relation, we 
observe that S±{uj)duj is the quadrature variance of the cavity signal-signal modes in the interval 






between uj and oj + dw. The quadrature variance in the interval between w' = — A and uj' = \ can 
then be expressed as 

(Aa±)L = y^5±(wW, (89) 

in which uj' = uj — ujq. 

We now proceed to obtain the spectrum of quadrature fluctuations for the cavity signal-signal modes 
produced by the degenerate subharmonic generator. In view of (1691) . one can rewrite Eq. ()87l) as 

Dp roo 

S±{oj) = — / (a±(t)a±(t + (90) 

TT Jo 

Moreover, upon setting ki = K 2 = k and applying in Eqs. (fT3]) and (fTTl) . we get 

^{di) =-^K{di) - £{al) (91) 

and 

^( 02 ) = -^^( 02 ) - e(al), (92) 

in which e is defined by (l30|) . Now addition of Eqs. (IMIl and (1921) results in 

-£{d^), (93) 

where d is defined by (l36l) . Hence using Eq. (1931) and its complex conjugate, one readily hnds 

^(a±(t)) = -r]±{d±{t)), (94) 

in which (a+(t)) and (a_(t)) are given by (l67)l and (l68l) and 

V± = ^±e. (95) 

We note that the solution of Eq. (I94jl can be written as 

{d±{t + t)) = {d±{t))e~'^^'^ (96) 

and applying the quantum regression theorem, we have 

(a±(t)a±(t + T)) = (a|(t))e-^±E (97) 


Therefore, on account of (f^ . one can put Eq. ([90]) in the form 


1 

S±iuj) = -{Aa±)^Re / (98) 

^ Jo 

so that on carrying out the integration, the spectrum of the plus quadrature fluctuations for the 
cavity signal-signal modes is found to be 


5 '+(a;) = (Aa+)" 


(f +£)A 


(w — ujoY + [f + 


and on taking into account d74)l . we get 


S+{u:) = 


(f + e)/ 7 r 


(a;-u;o)2 + [f+ e]2 


2k 


(99) 


K T 26 


( 100 ) 






Fig. 6.2 A plot of local quadrature squeezing versus A [Eq. (11051) 1 for e = 0.4 and k = 0.8. 


Furthermore, upon integrating Eq. (jlOOp in the interval between oj' = —A and oj' = A, applying the 
fact that 


fA 

-A 


duj' 


I -X 


-II A 
= —tan — 


we arrive at 


(Ao+)^;^ = —tan 


-1 


TT 


f+e 


2k 


K “1“ 2s 


On the other hand, on setting e = 0 in (I1U2I) . we find 


i^a+)l±\ = -tan 

TT 


_i/2AA 


. • 


( 101 ) 

( 102 ) 

(103) 


This represents the quadrature variance of the two-mode cavity vacuum state in the same frequency 
interval. We define the quadrature squeezing of the cavity signal-signal modes in the interval between 
uj' = —A and w' = A by 


5±a = 


(Aa+)^^;^ - (Aa+) 


±A 


Therefore on account of (|102p and (llOSp . we see that 


5±a = 1 - 


tan i(^) 


2tan \n + 2e) 


2k 


(104) 


(105) 


The plot in Fig. 6.2 shows that the maximum local quadrature squeezing is 74.9% below the 
vacuum-state level. This occurs in the frequency interval X± = 0.05. In addition, we note that the 
local quadrature squeezing approaches the global quadrature squeezing as A increases. 


6 Conclusion 

The mean photon number of the twin light beams with the same frequency, obtained applying the 
conventional Hamiltonian, is found to be half that of the twin light beams with different frequen¬ 
cies. This unexpected result must be due to representation of the twin light beams with the same 


















frequency by second-order annihilation and creation operators in the conventional Hamiltonian. We 
then realize that the twin light beams with the same or different frequencies must be represented 
in the Hamiltonian by first-order annihilation and creation operators. We thus describe the process 
of subharmonic generation leading to the creation of twin light beams with the same or different 
frequencies by a Hamiltonian of the form indicated in Eq. (HD- 

Evidently, the mean photon number calculated using this Hamiltonian turns out to be the same for 
both the signal-signal and the signal-idler modes. In addition, the mean photon number of the pump 
mode is found to be the mean photon number of the pump mode in the absence of the subharmonic 
generation minus the mean photon number of either the signal-signal or the signal-idler modes. 
This is exactly what we would expect the mean photon number of the pump mode to be. On the 
other hand, we have established that the maximum global quadrature squeezing of the signal-signal 
or the signal-idler modes is 50% below the vacuum-state level and the maximum local quadrature 
squeezing of the signal-signal modes is 74.9% below the same level. 


Appendix: The Conventional Hamiltonian 


Here we wish to calculate, employing the conventional Hamiltonian, the mean photon number of 
the twin light beams with the same frequency produced by a degenerate subharmonic generator. 
The process of subharmonic generation leading to the creation of the twin light beams is usually 
described by the Hamiltonian 

H = — b) + — 6a^^), (Al) 


in which a is the annihilation operator for the twin light beams, b is the annihilation operator for 
pump mode, A is the coupling constant, and /r is proportional to the amplitude of the coherent 
light deriving the pump mode. We may refer to a Hamiltonian of the form described by (Al) as 
second-order Hamiltonian. We consider here the case in which the twin light beams and the pump 
mode are in a cavity coupled to a vacuum reservoir via a single port-mirror. The equations of 
evolution of (a^a) and (a^) can be written as 

d 1 

— (a^a) = —i([a^a, H]) H— KTr[{2apa^ — a)ap — pa^a)a^a] (A2) 

dt 2 

and 

— (a^) = H]) + -KTr[{2apa^ — a^ap — pa^a)a^], (A3) 

tlu ^ 

in which k is the cavity damping constant. Now taking into account Eq. (Al) and the fact that 


and 


we easily get 


-KTr[{2apa} — a)ap — pa}a)a)a] = —K{a)a) 


Tr[{2apa^ — a^ap — pa^a)S?] = —K{a^), 


dt 


{a)d) = —K{a)d) — A(6^a^) — A(6a^^) 


(A4) 


(A5) 


(A6) 


and 


A(a2) = -K{a^) - 2\{ba)a) - X{b). (A7) 

The steady-state solutions of these equations are 

(atfl) = --{Pa^) --{ba^^) (A8) 

K K 

and 

2A - A - 

(a^) =-(6a^a)- (b). (A9) 

K K 

Now substitution of (24) into Eqs. (A8) and (A9) yields 

{d)a) = --{a^) - (AlO) 

K K 

and 

(a^) = - — {a^a) - (All) 

K K 

where e is defined by Eq. (30). Einally, on account of Eqs. (AlO) and (All), we find the mean 
photon number of the twin light beams to be 
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